Directed transport of Brownian particles in a double symmetric potential 
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We investigate the dynamics of Brownian particles in internal state- dependent symmetric and 
periodic potentials. Although no space or time symmetry of the Hamiltonian is broken, we show that 
directed transport can appear. We demonstrate that the directed motion is induced by breaking the 
symmetry of the transition rates between the potentials when these are spatially shifted. Finally, we 
discuss the possibility of realizing our model in a system of cold particles trapped in optical lattices. 
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In usual Brownian motion , thermal fluctuations act 
as a reservoir of energy which distillation within the sys- 
tem of particles results in random trajectories and molec- 
ular chaos on a microscopic scale and in symmetric spa- 
tial diffusion on a macroscopic scale. The problem of con- 
verting the energy of noise into deterministic mechanical 
work (i.e. rectification of the energy of noise) to induce 
directed motion was first addressed by Schmoluchowski 
and Brillouin and has attracted much attention dur- 
ing the past years (see Ref. (3] and references therein). 
This question is crucial in particular for understanding 
how molecular pumps work in biological environnements 
M or for designing new synthetic pumps for chemistry 
[5j ■ It is by now clear ^] that the necessary and generally 
sufficient conditions for noise rectification are : (i) break- 
ing the spatio-temporal symmetry of the Hamiltonian (to 
match the Curie principle 0) and (ii) non thermody- 
namic equilibrium situations (to avoid detailled balance 
and resulting zero current |3|)- 

Most studies rely on Brownian particles experiencing 
a spatially periodic potential u{x) plus eventually a time 
T- periodic force f{t). After Curie's principle 7], if both 
symmetry conditions 



u{x) — u{2xq 

f{t + T/2) = - 



-f{t) 



(1) 

(2) 



are fulfilled (xq being any given position), then no di- 
rected motion is possible. On the contrary, breaking one 
of these symmetries is generally sufficient to ensure a net 
motion of the particles J9j . Both cases of spatial sym- 
metry (P) breaking U and temporal symmetry Q 
breaking have been widely studied in a large vari- 
ety of systems and evidences for directed transport have 
been given. In these systems, the non-equilibrium re- 
quirement is ensured by either non-gaussian noise, tem- 
perature oscillations or flashing potentials In re- 
cent years, more elaborated systems that mimic biologi- 
cal system more properly have been investigated. It has 



been thus demonstrated that interestingly directed mo- 
tion can be obtained in symmetric potentials and in the 
absence of a time-dependant force, i.e. in systems where 
conditions CJ-© are fulfilled. This requires complex 
situations, e.g . position-dependent mobility in flashing 
potentials jfj, state-dependent mob ility jl2 ] or time- 
dependant spatial fluctations of walls . 

In this work, we propose a very simple system for which 
directed motion is obtained despite validity of both spa- 
tial and temporal (0) conditions. Here, it is the whole 
dynamics made of Hamiltonian motion phases and dis- 
sipative processes which is non-symmetric. More pre- 
cisely, directed motion is induced for low-damped parti- 
cles in internal state-dependent periodic potentials with 
spatially shifted pinning points and for different transi- 
tion rates between the potential curves. The behavior 
is explained on symmetry grounds and with the help of 
a simple rectification process that qualitatively describes 
the dynamics of the particles. Our conclusions are sup- 
ported by a study of the influence of the relevant param- 
eters (shift between the potentials and strengths of the 
transtion rates). We finally discuss the possibility of ex- 
perimental observation of such an effect in a system of 
cold particles trapped in optical lattices. 
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FIG. 1: Diagram of the system considered in this work. A 
Brownian particle can live in two different internal states 
and 1 corresponding to (^-shifted potentials; the transition 
rates are asymmetric (7+ 7^ 7„). The rectification process is 
shown in the right-hand side of the figure (see text below). 



'Present adress:Institut fiir Theoretische Physik, Universitat Han- 
nover, D-30167 Hannover, Germany 



Consider a Brownian particle living in N internal 
states (indexed hy j G [0..N — 1]) corresponding (eventu- 



2 



ally) to different momentum diffusion constants D^j{x) 
and different spatially modulated potentials Uj{x) and 
submitted to a time-dependent homogeneous force fj {t) . 
For example, the different potentials may correspond to 
different chemical species and the transitions to chemical 
reactions or to different quantum internal atomic states 
and light-induced transitions respectively. 

The Fokker-Planck equation (FPE) for the population 
Wj{x,v,t) of the state characterized by x,v (external 
state) and j (internal state), reads 

[dt + vd,] W, +d,[v + u'^{x) - fj{t) + D,j{x)d,] W, 

= E ^r^j {^)Wr - E 7.-/ {x)W, (3) 

where ^j'^j represents the transition rate from internal 
state j' to state j. The FPE Q is written in time units 
of the friction coefficient and in space units of the typical 
variation scale of the potentials Uj so that all variables 
are dimensionless. The left-hand side of equation Q is 
a usual FPE for a particle submitted to a friction force, 
a potential-induced force, a time-dependent force and a 
momentum diffusion and thus captures the one dimen- 
sional Brownian motion in a regime ranging from over- 
damping to low damping. The left-hand side of equa- 
tion is a source term that accounts for random jumps 
between the internal states. Non-thermal transition rates 
{e.g. induced by chemical reactions or optical transitions) 
result in non-thermodynamic equilibrium situations 'sj. 
In the following, we suppose that 

- Uj is spatially periodic and symmetric in the sense 
of O, 

- fj is T- periodic and fulfills lO; note that the aver- 
age value of fj vanishes: i JJ" dtfj{t) = 0. 

We consider the case of a double symmetric potential 
{N = 2) with ((/)-) shifted maxima (see Fig. ^| with (in 
general) non-symmetric transition rates (70,1 7^ 7i,o) a-nd 
in the absence of a time-dependent force {fj = 0). More 
precisely, we use sine potentials with the same period 27r: 

U(){x) = cos(a;) 

ui{x) = AiCoa{x + (f>) . (4) 

We assume that the mean kinetic energy (proportion- 
nal to the temperature) of the system is smaller than 
the depth of the trapping wells corresponding to the po- 
tential minima {i.e. typically < Aj). Moreover, we 
assume the energy damping to be much smaller than 
the typical frequency of the oscillations in the potential 
wells {i.e. y^\Aj \ 3> 1) so that the dynamics falls into the 
low damping regime. Thus the particles are essentially 
trapped in the wells and their basic behavior consists in 
weakly damped and noisy oscillations at the bottom of 
the potential minima. 

For the sake of simplicity, we finally assume minimal 
requirements for the phenomenon of directed transport 
demonstrated in the reminder of the paper: 



- Ao = Ai = A, 

- Dvj is position- and internal state- independent 
{D^, = Dv), 

- 7+(= 70^1) a-nd 7-(= 71^0) are position indepen- 
dent and different in general. 

Focus first on the Hamiltonian part of the dynam- 
ics. It is clear that the Hamiltonian is symmetric with 
respect to both space and time. Indeed, it is invari- 
ant through {ue{x)} {ui-e{2xo — x)} where xq is the 
middle point between a maximum of uq and a maxi- 
mum of Ml. Moreover, there is no explicit time de- 
pendence in the equations of motion Therefore, 
the Hamiltonian part of the dynamics is clearly sym- 
metric with respect to both conditions ^ and ^ and 
cannot induce any directed motion. However, the dy- 
namics has to be regarded as an intercombination of 
Hamiltonian motion (motion on the potential curves) 
and dissipative processes (transitions between different 
internal states). Here, the system no longer preserves 
the spatio-temporal symmetry since it is not invariant 
through {uc(x), 7e} ^ {wi-c(2a;o — a;), 7i_c} except for 
74. = 7- . Hence different transition rates break the 
spatio-temporal symmetry of the system and force-free 
directed motion can be expected. It is the aim of the fol- 
lowing paragraphs (i) to demonstrate this phenomenon 
and (ii) to describe the corresponding microscopic mech- 
anism. 

We performed numerical simulations of the dynamics 
of the particles in the system described above by means 
of a Monte-Carlo integration with a set of typically 2000 
atoms. It is then straightforward to determine the mean 
position of the cloud of particles as a function of time 
from which we deduce the average velocity {V) of the 
cloud: 



where ( . ) stands for the ensemble average over all parti- 
cles. For all the results presented in this work, we found 
a constant mean velocity (see inset of Fig.|2l. 

We first study the behavior of the mean velocity V as 
a function of the space shift (j) between the potentials uq 
and ui in a situation where the transition rates between 
the internal states are significantly different (7- = 87+). 
The results of the numerical simulations are shown in 
Fig. [2 The cloud of particles clearly shows a directed 
motion through the potentials for almost all values of 
(j). These results show that a directed transport of low- 
damped Brownian particles can be induced in a 2 state 
spatially symmetric potentials with shifted minima. The 
main features of Fig. |21 can be explained under sym- 
metry considerations. First, the average velocity van- 
ishes {V = 0) only for discrete values of (f> correspond- 
ing to symmetric situations {(j) = 0, tt). Here the symme- 
try is different from the one discussed above and reads 
{u^{x),j^} <-> {ue(— x),7e}. Notc that this space inver- 
sion is valid only for (j) = or (j) = tt and thus V = only 
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for these values of (j>. Second, changing (p into {2tt) — 0, 
the shift of the minima of ui with respect to the minima 
of uq gets opposed and this results in the inversion of 
the direction of motion of the particles as observed on 
Fig. 12 Third, we find that the mean speed {\V\) of the 
directed transport is maximum for values of close to 
7r/2 or 37r/2, i.e. in a situation where the minima of the 
potentials are almost maximally shifted. 




FIG. 2: Mean velocity of the particles as a function of the 
phase 0. Inset: displacement of the particles as a function 
of time for two different values of (p. Numerical parame- 
ters for this simulation are A = 200, Dv = 75, 7- = 0.75r, 
7+ = 0.25r and F = 10. 

We now describe the noise rectification mechanism. In 
the low damping regime considered in this work, the 
microscopic dynamics of the particles consists in noisy 
oscillations at the bottom of the potential wells inter- 
rupted by random transitions from a potential curve to 
another. Suppose now that < 7+ < 7-. Then a par- 
ticle spends the most of time in potential uq and oscil- 
lates around a given minima. Due to the finite value 
of 7-1- , it may be transfered into potential ui. If the 
minima and maxima of uq and ui are spatially shifted 
{(j) 7^ riTT, with n an integer), the particle initially local- 
ized near a minimum of uq falls away from a minimum of 
Ml and consequently experiences a force oriented towards 
the left (for 2mT < < (2n + Ijir) or towards the right 
(for (2n — l)7r < < 2n7r). It thus moves in the corre- 
sponding direction. It is finally rapidly pumped back to 
potential uq, generally in the initial trapping well. As a 
consequence of this cycle (see right hand side of Fig.^l, 
a mean directed force is exerted on the particle and the 
center of the oscillations of this particle in a well of uq is 
shifted (towards the left for 2n7r < (j> < {2n + l)7r or to- 
wards the right for (2n — l)7r < < 2n7r). This formally 
lowers the potential barrier for escaping the trapping well 
in the corresponding direction and rises it in the opposite 
direction. This clearly induces a directed average mo- 
tion of the cloud of particles. Note that this mechanism 
matches the mean features of Fig. |5| described above. 

As discussed previously, considering different values 
of the transition rates 7+ 7^ 7_ is crucial to break the 
spatio-temporal symmetry of the system and correspond- 
ingly to induce directed motion of the particles through 
the periodic structure. We now investigate the infiuence 



of the transition rates on the average velocity V. We 
plot y as a function of 7_|_ on Fig. O for three differ- 
ent values of the total transition rate F = 7+ + 7_ and 
in a situation where the potentials uq and ui are sig- 
nificantly shifted {</) ~ 37r/2). Wc find that V strongly 
depends on 7+/F as expected. For 7^/F = (7-/F ~ 1) 
or 7_|_/F = 1 (7-/F = 0), the system is equivalent to a 
Brownian particle in a single symmetric potential and 
therefore does not show any directed motion: F = 0]- 
7+/F = 0.5 (7+ = 7_) corresponds to a symmetric situa- 
tion as discussed above and results in no directed motion 
as confirmed by Fig.O For all other values of 7±/F, a di- 
rected motion is clearly demonstrated. The maximum of 
y| is found to correspond to strongly asymmetric values 
of 7-1- and 7_ (typically for 7±/7ip — 8). This is consis- 
tent with the rectification mechanism proposed above. It 
is indeed clearly necessary that particles get trapped in 
the wells of (e = 0, 1) and encounter short transitions 
in wi-e to induce significant directed motion. We finally 
note that 7+/F 1 — 7+/F results in V —V. This is 
an obvious consequence of the anti-symmetry of the sys- 
tem through exchange of the internal states 0^1 and 
corresponding change of the phase shift: (/)<->—(/). We 
finally note that \V\ decreases when F increases and this 
indicates that the rectification process is no longer effi- 
cient for high values of F. Indeed, in this case, a particle 
undergoes many cycles so that it experiences an average 
potential (7-Uo + 7+ui)/F which is symmetric and again 
no directed motion is possible |2j|- 




FIG. 3: Average velocity V of the cloud of particles as a 
function of the transition rates between internal states. V is 
plotted versus 7+ for a fixed value of F = 7+ -f- 7_ . Parame- 
ters of the simulation are A = 200, Dv = 75 and (j) = 3n/2. 

We now discuss the possibility of realizing our model 
in experimental setups. Ratchet effects have been stud- 
ied in a large variety of systems 3] but few of them al- 
low for independant control of the physical parameters in 
wide ranges with accuracies as good as in systems of cold 
atoms trapped in optical lattices 15, 16]. We thus pro- 
pose to use such a system to observe the phenomenon of 
directed transport in a double symmetric potential. Opti- 
cal lattices are formed by the interference pattern of laser 
beams In appropriate configurations, matter- light 

interactions create (i) a periodic and symmetric potential 
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for the atoms, (ii) a friction-like force and (iii) random 
fluctuating forces. The competition between the last two 
phenomena results in an equilibrium kinetic energy lower 
than the potential depth. The atomic dynamics has been 
shown to be equivalent to Brownian motion in a periodic 
potential [3) ^^^'^ low-damping regime [23| . There- 

fore these properties match the requirements considered 
in this work for observing directed transport in a double 
symmetric potential. In fact, the parameters used for the 
simulations presented in this work roughly mimic optical 
lattices. In order to induce a double potential, we sug- 
gest to follow the scheme recently presented in Ref [21| . 
By illuminating ^'^^Cs with two lasers respectivelly near- 
resonant with two atomic transitions, the atoms see a 
double sine-like potential. Moreover, the atom-light in- 
teraction induces non-thermal transition rates between 
these potentials. It is finally worth noticing that the 
space shift between the potentials, the potential depths 
and the transition rates can be controlled experimentally. 
Therefore, we believe that this system matches the main 
characteristics of the model studied in this work (phase- 
shifted symmetric potentials, non-thermal and different 
transition rates, temperature below the potential depth, 
low-damping regime) and should be well suited for ob- 
serving directed motion. 

To sum up, we have studied the directed motion of 
particles induced by spatio-temporal symmetry breaking 
in a multi-level periodic system. Our system consists in 
low-damped Brownian particles in a multi-level flashing 
periodic and symmetric potential. We have shown that 
directed motion is induced whenever the minima of the 



different potentials are shifted and for different transi- 
tion rates. This breaks the spatio-temporal symmetry of 
the complete dynamics of the particles which corresponds 
to a combination of Hamiltonian motion (on the poten- 
tial curves) and dissipation (related to the internal state 
transitions). We have finally proposed a system based on 
cold atoms trapped in an optical lattice where directed 
transport in a double symmetric potential may be exper- 
imentally observed and characterized with a high degree 
of accuracy. 

Exploring Brownian motors and corresponding ratchet 
effect helps in understanding the physics of protein mo- 
tors, Josephson junction arrays and quantum dots. Usual 
models consider particles in space asymmetric potentials 
or undegoing time-asymmetric forces More re- 

cent models consider elaborated schemes in symmetric 
potentials |TlL IT^ fl^. In this work, we have proposed 
a very simple scheme in which ratchet effect can be ob- 
served in symmetric potentials. In fact our system con- 
siders a low-damped Brownian particle minimally in a 
2-state potential with shifted minima and different tran- 
sition rates. This simplicity makes the scheme quite gen- 
eral and it could be observed in various systems for ex- 
ample in cold atoms optical lattices. 

We are indebted to Professor Anders Kastberg for en- 
lightening discussions on the physics of double optical 
latices and for valuable comments on the paper. Labora- 
toire Kastler Brossel is an "unite mixte de recherche de 
I'Ecole Normale Superieure et de I'Universite Pierre et 
Marie Curie associee au Centre National de la Recherche 
Scientifique (CNRS)" . 
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